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Abstract 

An analysis is presented of the suitability of the method of equivalent 
linearization for estimating the fatigue life of a nonlinear structure. Comparisons are 
performed of the fatigue life of a nonlinear plate as predicted using conventional 
equivalent linearization and three other more accurate methods. The excitation of 
the plate is assumed to be Gaussian white noise and the plate response is modeled 
using a single resonant mode. The methods used for comparison consist of 
numerical simulation, a probabilistic formulation, and a modification of equivalent 
linearization which avoids the usual assumption that the response process is 
Gaussian. Remarkably close agreement is obtained between all four methods, even 
for cases where the response is significantly nonlinear. 

Introduction 

Equivalent linearization has been applied to predict the mean square 
responses of a variety of complex nonlinear systems. However, the applicability of 
the method for the prediction of other response parameters such as the power 
spectral density and higher order statistics is severely limited by the assumption that 
the probability density of the system response is Gaussian. It is known that the 
response of a nonlinear system with Gaussian input will not be Gaussian. It has 
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been shown recently 1 that the response spectrum predicted by equivalent 
linearization for a Duffing oscillator differs substantially from that obtained by 
numerical simulation. Numerical simulations have shown that as the system is 
excited by higher level inputs, the resonant response peak, observed in the power 
spectrum, tends to broaden. Conventional equivalent linearization does not predict 
this effect. A modification to equivalent linearization has been developed, however, 
which avoids the assumption that the response is Gaussian and accurately predicts 
the broadened spectral response at high excitation levels . 1 Since conventional 
equivalent linearization is not suitable for predicting spectral response, there is a 
concern that it may not be useful for fatigue life estimation. In the study described 
herein, it is shown that, in spite of its shortcomings, conventional equivalent 
linearization provides fatigue life estimates which agree closely with those obtained 
by more sophisticated methods for a system which vibrates nonlinearly in a single 
resonant mode. 

The estimation of fatigue life typically depends on an accurate description of the 
response probability density and the response spectrum. Both of these parameters 
influence the probability density of response peaks which is a key ingredient in 
fatigue estimation. Both of these factors are very crudely represented in 
conventional equivalent linearization. The method remains quite attractive, however, 
since a detailed analysis of the fatigue life of a complex nonlinear structure using 
more accurate methods such as numerical simulation or probabilistic techniques is 
not presently feasible unless the structure has a very limited number of degrees of 
freedom. If equivalent linearization could be applied directly to estimate the fatigue 
life of a complex nonlinear structure, it would avoid either the extensive numerical 
complexities of time domain simulations, or the mathematical intricacies of obtaining 
a detailed probabilistic description. 
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In the following, an overview will be presented of basic fatigue life estimation for 
a structure which experiences random loads. 'Classical' analytical methods will then 
be discussed which may be used to obtain fatigue life estimates for linear and 
nonlinear structures having a single degree of freedom. A numerical simulation 
approach is then presented. The application of conventional equivalent linearization 
for estimating fatigue life is then described. Finally, estimates of the fatigue life are 
presented using a modified equivalent linearization approach which avoids the 
Gaussian assumptions. 1 These techniques are then applied to predict the fatigue 
life of a nonlinear plate driven by Gaussian white noise. Remarkably good 
agreement is obtained between all methods. These results indicate that it may be 
possible to adapt equivalent linearization to predict the fatigue life of complex, multi- 
mode, nonlinear systems. 

Fatigue Life Estimates for a System with Random Stresses 

For the present study we will assume that the relation between the stress 
amplitude, S, and the number of cycles to failure N, may be approximated by the well 
known relation, 



where c and b are experimentally obtained constants for a given material. It will be 
assumed that the number of fatigue cycles experienced by the structure is equal to 
the number of positive stress peaks, or stress reversals, that occur overtime. A 
sample time history is shown in Figure 1 where the positive peaks are identified by 
circles. 

The assumption that damage occurs at each positive stress peak may be 
conveniently combined with the Palmgren-Miner linear damage accumulation rule 2 
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to produce fatigue life estimates. In this theory, the damage, Dj, caused by stress 
reversals at the stress level Sj is 


n(S.) 
D i (S i )= N(S.) 


( 2 ) 


where n(Sj) is the number of stress reversals experienced by the structure at the 
stress level Si and N(Si) is the number of reversals required to cause failure at this 
stress level. The total damage, D m , will be the sum of the damage at all stress levels 
that occur, 


n(SJ 


D m =SD i (S i )=I N(ST 
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( 3 ) 


From equation (1), 
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( 4 ) 


so that equation (3) becomes 

D m =4 s "( s i> s i 


( 5 ) 


The form of equation (1) assumes that all stress peaks occur at positive stress 
levels as would be the case for sinusoidal loading. When a resonant system is 
subjected to random loading, however, we must allow for the possibility of stress 
peaks at negative stresses which make positive contributions to the accumulated 
damage. To account for this, we will take the absolute value of the stress level, Si, so 
that equation (5) becomes 



7S n(S ) 


S. 


( 6 ) 


Failure is predicted to occur when D m =1. 

Since the response and rate of damage accumulation in the structure is 
assumed to be random, the fatigue life may be estimated from the expected value of 
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the rate of damage accumulation. This may be expressed as an integration over all 
possible stress levels in the form 


oo p 

E[D m |= J n(S) J|-P(S) dS 


(7) 


where P(S) is the probability density for response peaks. 


The assumption that the response process is stationary gives the fatigue life as 

t= TO ' (8) 

In the case where the system is linear and has a single lightly damped resonant 
mode, the peak probability density may be approximated by a Rayleigh distribution 
for a narrowband process, 


P(S) = ^ 
o „ 


2a ‘ 


S>0 


(9) 


where ^ is the mean square stress level. 

Approximate Methods for Single Degree of Freedom Systems 

For a linear single degree of freedom system, Miles 3 has replaced n(S) in 
equation (7) with the system resonant frequency in Hertz, co 0 /27t. Equations (7) and 


(9) then give Miles' single degree of freedom formula 3 , 

E [ D m] = i^ (V^o s ) b r (b + z 


where, 


Hy) = 2 J 


2y-1 -x . 
x e dx, y > 0 


( 10 ) 


( 11 ) 


If a single degree of freedom oscillator is nonlinear it is possible to evaluate 
equation (7) by replacing n(S) with the expected total number of peaks per unit time, 
E[Mj]. If we again assume that the response process is narrowband, it is not difficult 
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to approximate the peak probability density, P(S). If the displacement response, x(t) 

is governed by Duffing's equation, 

x+ cOq (x + ex 3 ) + ax = f (t ) , ^ 

where co 0 is the linear natural frequency, e and a are constants, and f(t) is Gaussian 
white noise, then the damage rate may be approximated by 4 , 

2 4 

X EX 

00 2 2 

1 o b f h ^ 2 

E [ D m]) = c C 0) o K J x °( x + ex 3 )e 0 dx , (13) 

0 


where 


1 




2 2 
™ 0°0 




dx 


(14) 


and o^is the mean square displacement response when e=0, and 

= Gf/(4<Uoa), where Gf is the constant single-sided power spectral density of 

the excitation in units of f 2 /Hertz. The constant K in equation (13) is an assumed (for 
the present study) linear relation between stress and displacement, i.e. 

S = K x (15) 

When e=0 equation (13) reduces to equation (10). 

An alternative approach to estimate the fatigue life is to solve equation (12) 
numerically and use equation (15) to determine the stress as a function of time. 
When a stress reversal occurs at the stress level S, one simply increments the 
accumulated damage according to 

D = D + Is£ 

u m c ’ (16) 

where D m is initially set to zero. If the solution is simulated for a time period equal to 
t then the damage rate is D m / * and the predicted fatigue life is 



( 17 ) 
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While the above methods are straight-forward to use for simple systems with a 
single degree of freedom such as those governed by equation (12), there are 
substantial difficulties in applying them to more complex multi-degree of freedom 
structures. The derivations of equations (10) or (13) rely on the assumption that the 
response process has a very narrow bandwidth, which is expected to be true only for 
a lightly damped single degree of freedom oscillator. Numerical methods such as 
equation (16) are computationally impractical when simulating the random response 
of a complex nonlinear system. Since most structures of practical interest are too 
complicated to analyze using existing techniques there is a need to develop more 
powerful procedures. 

Fatigue Life Estimates Using Equivalent Linearization 

In this section, methods for estimating fatigue life based on conventional 
equivalent linearization 5 and by a modified approach will be described. The 
modified equivalent linearization method discussed here follows from reference [1] 
where it is shown that this technique gives very accurate estimates of the power 
spectral density of the response of a nonlinear system which is a dramatic 
improvement over the conventional approach. It will be shown that both of these 
theories may be adapted to give reasonable fatigue life estimates for a system with 
nonlinear stiffness that has a single resonant mode. 

The method of equivalent linearization consists of approximating equation (12) 
by 

x + cOgX + cxx + e(x) = f(t) , 

where co e is the equivalent linear natural frequency and e(x) is the error in the 
approximation. co e is chosen so that e(x) is minimized. The main difference between 
the conventional approach and the modified method 1 is that in the modified method, 
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03 e is taken to be random rather than deterministic. One may then avoid the 
assumption that the response is a Gaussian random process. 

In conventional equivalent linearization 5 the square of the equivalent linear 
natural frequency is given by 

“ 2 e = lf( 1 + V 1 + 12 “o) ' 09) 


where, as before 



4(£> 2 Q a 


( 20 ) 


is the mean square response when c=0. 

The mean square value of the response of the nonlinear system is then 
approximated by 


4co„a 


( 21 ) 


To obtain an estimate of the fatigue life of a system governed by equation (12) 
we could simply calculate the fatigue life of an equivalent linear system such as in 
equation (18) with e(x) neglected. This may be accomplished by using equation (10) 


with co 0 replaced by co e and o s replaced by K o x , 


to b / 

E [ D -.]=^ K «.) r \ 


b + 2 


( 22 ) 


This result is derived with the assumption that the peak probability function is a 
Rayleigh distribution which corresponds to a narrow-band Gaussian random 
process. This assumption may be avoided when the modified method of reference 
[1] is used. In reference [1] the system is considered to be excited by a driver having 
a randomly varying power spectral density. Since the excitation is Gaussian, the 
input spectrum level may be described in terms of a chi-square random variable, k. 
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The random fluctuations in to 2 are then described in terms of where £ is 
exponentially distributed, P(^)=e - ^, 4>0, and 



where k is obtained by iteratively solving 



In reference [1] expected values are calculated in terms of integrations over 
P( 4 ) rather than in terms of the response variable, x, as in the conventional 
formulation. Since the probability density of x is usually unknown, a Gaussian 
distribution is typically assumed. Equation (7) may be transformed into an 
integration over as 

E [°m] = J n <%> ^ p tt> • (25) 

o 

From equation (19) we may set 



In reference [1] a result is presented for the 'mean-square' response, c 2 (£), 


corresponding to each value of 



(27) 


In equation (25) we require an expression for the peak stress as a function of % . An 


approximation for the peak response level, IX p (£)l may be obtained from equation 
(27) by considering the response corresponding to a particular value of % to be 
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dominated by a single frequency component at the frequency o) e (^). With this 
assumption the relation between the mean square and the peak response is 




,Kf 


or |x p (E.)|= V 2 °. 2 © 


From equation (15) and (28) the peak stress m ay be written as 

|SK)| = K <%&) • 


( 28 ) 

( 29 ) 


Equations (25) through (29) then give 

b-1 3zb . 

2 VkV T f I ~ \ 2 2 -t 

E t D m] = 2 rc H 1 + V 1+12Ek °oU 5 e ' 5 ( 30 ) 

0 

The fatigue life is obtained from equations (8) and (30). When e=0 , equation (30) 
reduces to equation (10) for the linear system. 

Governing Equation for a Nonlinear Plate 
To compare estimates of the fatigue life of a system using the methods of the 
previous sections we will calculate the coefficients of Duffing's equation (12) 
corresponding to the response of a nonlinear plate with simply-supported 
boundaries. We will assume that the plate response may be described by a single 
resonant mode. Since we are only attempting to compare solution methods it is 
sufficient to use an approximation to obtain the coefficients of equation (12). It is 
important, however, to perform the comparisons in a region of parameter space that 
is physically relevant. 

If the in-plane motion of the plate is constrained at the edges, the Berger 
hypothesis 6 leads to 

4 2 

DV w-ph e C p (w xx + w yy ) + phw + aw = p(x,y,t) (31) 


where 


*" y ^ x 


e = 


1 


2L x L v 
x y 0 0 


/ J t" 2 x + w 2 y )dxdy 


(32) 


w is the transverse displacement, p is the density, h is the thickness, a is a viscous 
damping coefficient, and p(x,y,t) is the applied pressure. L x and L y are the 
dimensions of the plate in the x and y directions. The constants D and C p are the 


bending rigidity and wave speed, 


D = 


E h 




12(1 -if) r p( 1 — u / 

E is Young's modulus and u is Poisson's ratio. Comma denotes partial 

differentiation. For a simply supported plate, 

w(x,y,t)= A(t) sinf-p-^ sir|p^ 


(33) 


(34) 


V I) 


where A(t) is an unknown function of time. 

By substituting equation (34) into equations (31) and (32) the governing 


equation takes the form of equation (12), 

2 . - - 3 . 


A+co q (A + eA ) + aA = F(t) 


where, 


2 __ _D_ 
w o- ph 


K 


K 


e = 


2h 


(35) 

(36) 

(37) 


and 


L y L x 


1 1 s, i s K§) p(x,y,t) dxdy • (38) 


Let the excitation have the form 


P(x.y,t) = p 0 (t) sir^-p-j sin 


( 

, L y. 


(39) 


Equations (38) and (39) give 
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F(t) = 


P 0 (t) 

ph 


(40) 


The damping coefficient, a, in equation (35) is taken to be 

a = wo n , (41) 

where q is the loss factor. 

For the present study we will obtain the stress-displacement relation, equation 
(15), by neglecting in-plane displacements and assuming that the maximum stress 
may be approximated by 

S = E e y , (42) 

where E y is the strain in the y direction at the surface of the plate at its mid-point 
(x = L x /2, y = Ly/2). For a simply-supported plate this will be the location of 
maximum stress. It is assumed that the y component of the strain is greater than the 
x component which will be true if L x > L y . Evaluating the stress in this case gives 


S = - E 4 % 
2 


E A sin 2 ^ 


2 L 


Eh 


V ) 


( k \ 

< L I > 


(43) 


The constant K in equation (15) is then 

El 
4 


y _ Eh ( n 


v L y 


(44) 


This is, of course, a very crude approximation but it will suffice for our present 
purposes. 

The remaining constants to be determined for our comparisons are the material 
properties b and c in the damage accumulation model, equation (6). These must be 
obtained experimentally. For this study, we only need to use 'typical' numbers for a 
common material. The constants used in the calculations are 

b = 6.33 and c = 6.56 x 10 30 . (45) 

The results presented in the following correspond to an aluminum plate with the 
following properties: (English units) 
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E = 1 0 7 PSI, p = 0.1 Ibs./in 3 , \> = 0.3, L x = 20 in., 

L y = 10 in., h = 0.063 in., and rj = 0.05. 

The excitation is assumed to result from acoustic pressure fluctuations. The 

driving sound pressure level is given by 

L =20 log -jp- ( 46 ) 

K H ref 

where p re f = 0.0002 p bar. A pressure of one pound per square inch corresponds to 
a sound level of approximately 170 dB. 

Results and Recommendations 

The calculated results are given in figure 2 for a range of excitation levels from 
130 to 155 dB. The predicted fatigue lives based on damage rates estimated using 
the 'classical' result, equation (13), numerical simulation as in equation (16), and the 
methods based on equivalent linearization, equations (22) and (30), are found to be 
in remarkably close agreement. For comparison, the fatigue life obtained for a linear 
system (e=0) is also shown. Since the nonlinearity limits the response amplitude in 
this case, the linear system has a shorter fatigue life. This may not be true in general, 
however, since in the present study a linear stress displacement relation (equation 
(15)) has been assumed. The effects of nonlinear in-plane strains should be 
included in future studies. 

The accuracy of equivalent linearization methods in predicting the fatigue life 
when the response deviates significantly from the linear behavior as shown in figure 
2 is quite remarkable. The method of equivalent linearization provides the 
parameters of a linear system which has 'nearly' the same mean-sauare response 
as the nonlinear system being studied. However, higher order statistics of the output 
of the equivalent linear system which are thought to strongly influence fatigue life 
(such as distributions of peaks, level crossing and spectra) differ substantially from 
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those of the nonlinear system. In this case, as often occurs in nonlinear analyses, 
the range of validity of the approximate solution methods are far from obvious and 
can be quite difficult to determine. The numerical results presented here indicate 
that the fatigue life of a nonlinear system with a single degree of freedom depends 
primarily on the mean square response level and the average frequency. A more 
detailed analysis is needed to verify this result. 

An examination of the error in the fatigue life as predicted using equivalent 
linearization could proceed by comparing the equations for the expected value of the 
damage rate, equations (19) through (22) with the probabilistic formulation in 
equations (13) and (14). It should be possible to obtain an expression for the 
difference between equations (13) and (22) which would explicitly show how the 
parameters of the system influence of the error. This analysis might lead to an 
improved approximate formulation in which the fatigue life is estimated directly rather 
than needing to first linearize the system. An approach of this type would be 
beneficial since it would be easier to see how the fatigue life is influenced by the 
system properties. 

While a goal of this investigation has been to improve our understanding of the 
fatigue of complex aerospace structures, the present study has been limited to 
systems which have a single degree of freedom. However, measurements of the 
response of stiffened fuselage structures can exhibit a nearly uncountable number of 
resonant modes. This is illustrated in figure 3 which shows the power spectral 
density of the strain at the location of maximum strain in a realistic test structure. 7 
The structure consists of a panel with orthogonal stiffeners such that there are a total 
of nine skin panels. The figure shows a large number of spikes which presumably 
correspond to the response of resonant modes. In acoustic fatigue studies it is 
common practice to treat structures such as that in figure 3 as a single degree of 
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freedom system. In reference [7] this structure was considered to have only one 
resonant mode with a frequency of 380 Hz. It may be argued that this approach 
works because it has lead to successful designs. It is, however, far from clear (at 
least to the author) why this works. Considering the structure to be a single degree 
of freedom system, which it obviously is not, is not at all appealing and could yield 
incorrect designs. In order to develop reliable acoustic fatigue design procedures for 
future materials and structural configurations there is a need to extend the existing 
analysis methods to include many degrees of freedom. From the results presented 
here, equivalent linearization appears to be a promising method which could be 
adapted for fatigue life predictions of nonlinear structures with multiple modes. The 
present numerical study and the suggested analysis discussed above for the single 
mode case should be repeated with multiple modes included. 

Conclusions 

The results presented here show that equivalent linearization produces 
estimates of the fatigue life of a nonlinear plate which agree closely with those 
produced by other methods that describe the statistics of the response much more 
accurately. The present analysis is limited, however, to single degree of freedom 
systems that may be described by Puffing's equation. Further analytical work is 
needed to examine the applicability of equivalent linearization in fatigue life 
prediction and in the development of approximate methods for predicting the fatigue 
of complex nonlinear structures with multiple degrees of freedom. 
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Figure 2. Predicted fatigue life of a nonlinear 
plate. 
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Figure 3. Strain spectrum for panel s. 
(From Reference (7)). 
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